#include "uhm.hxx"

#define N 1

using namespace uhm;

typedef Solver_Base_ Solver;

int main (int argc, char **argv) {

  initialize(argc, argv);

  // ** Environments
  Int_
    type = Mat_Base_::FLOAT,
    uplo = Mat_Base_::FULL,
    sym  = Mat_Base_::SYM,
    rhs  = 1,
    blk  = 256;

  Int_
    kind = Mesh_Base_::NODE_DEFAULT;

  // ** Solver initialization
  Solver solver(type, uplo, sym, rhs, blk);
  
  std::cout << solver ;

  solver.report(std::cout);

  finalize();

  return true;
}


